A PULSE FISHERY MODEL WITH CLOSURES AS FUNCTION 
OF THE CATCH: CONDITIONS FOR SUSTAINABILITY 



FERNANDO CORDOVA-LEPE, RODRIGO DEL VALLE, AND GONZALO ROBLEDO 



Abstract. We present a model of single species fishery which alternates closed 
seasons with pulse captures. The novelty is that the length of a closed season 
is determined by the stock size of the last capture. The process is described by 
a new type of impulsive differential equations recently introduced. The main 
result is a fishing effort threshold which determines either the sustainability of 
the fishery or the extinction of the resource. 



1. Introduction 

1.1. Preliminaries. This work proposes a model of manageraent of closed seasons 
(also named seasonal closures) for fisheries. The FAO's fisheries glossarjQ defines 
a closed season of a single marine species as the banning of fishing activity (in an 
area or of an entire fishery) for a few weeks or months, usually to protect juveniles 
or spawners. The closed seasons are implemented by a regulator authority in order 
to limit the productive factors in certain area during a specific time interval. 

A fishery process could envisage alternated periods of closures and open seasons, 
where the fishing is allowed. The closures are introduced according to bioeconomical 
necessities and its scope considers several spatio-temporal possibilities, leading to 
concepts as: seasonal closures-no area closure, short term area closure-no seasonal 
closure, short term time and fully-protected area, see e.g., [13] for details. 

This article will consider the special -but important- case of bioeconomic pro- 
cesses having the following property |31j : low frequency / short duration of open 
seasons combined with large magnitude of capture. In this framework, we propose a 
feedback regulatory policy which defines the lenght of the next closure as a function 
of the present captured biomass. A consequence of this regulation is that under 
certain fishing effort (this concept will be explained later) threshold, the conver- 
gence towards a constant length closures ensuring the ecological sustaintability of 
the resource is obtained. 

The literature shows many cases of a long-term time and area closures. An 
example is the pacific herring [Clupea pacificus) fishery, where the regulatory mea- 
sures have included very short open seasons as two hours joined with other inputs 
and outputs restrictions [23 • Another example is given by a village-base manage- 
ment program in Vanuatu, where the stocks should be harvested in a sequence of 
brief openings interspresed with several years closures (see [IB], [SI for details). 
Cases of annual fisheries combining short periods with large captures, followed by 
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Figure 1. Captures of anchovy and conimon sardine (1991-1995) 
in southern Chilean coast [11] . The period with the largest fishing 
mortality is the summer, while the lowest one corresponds to the 
closed season. 
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Figure 2. Idealized scheme of the previous figure: intense fishing 
season alternated by absolute closures of nine months. If the lenght 
of openings are small compared with closures and the harvested 
stock is large, it can be considered as a pulse. 

(comparatively) low ones at the rest, are registered. For instance, the fishery of 
common sardine (Strangomera bentincki) and anchovy (Engraulis ringens) in the 
southern Chilean coast has a behavior described by the Figure 1. 



The problem of determining the lenght of a closure (in a context of short-term 
open seasons with intense fishing effort) ensuring the bioeconomical sustantability 
is a complex issue: indeed, short term closures followed by intense captures could 
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induce an overexploitation of the resource. On the other hand, long-term closures 
could have some unexpected drawbacks as: i) Bio-sustainable economic rent with 
negative average [B]. ii) Promotion of negative indirect effects [5], |17| when the 
resource is a top predator. An example is given by the closure (1989-1992) of the 
nioUusk Concholepas concholepas in the Chilean coasts [27 ■ iii) The phenomenon 
known as race for fish: fishing units try to outdo one another in fishing power and 
efficiency during the brief openings [7,^. 

1.2. Mathematical modeling. In this paper, we will assume that the fishery 
process has two different time-scales. In the first scale (closed seasons), the growth 
of a single unstructured marine resource is described by using ordinary differential 
equations. Nevertheless, as the capture has a duration considerably shorter than the 
closures lenght, every open season will be considered as an instant of capture, i.e., 
a term of a sequence {tk}k, this is the second scale (See Figure 2). In consequence, 
the global process will be described by an impulsive differential equation, IDE (see 
e.g., [H], [24] for details). 

IDE equations have been used in the mathematical modeling of processes involv- 
ing impulsive harvesting: e.3.,[5],[l],[S],|Tl],[3g,[33],|31],E3 (see [22, and [3l] 
for more applications of IDE to bioeconomics and ecological processes respectively). 
In all these references, the following sequence of harvest instants is employed: 

(1.1) tk ~ kr with fc = 0, 1,..., and r > 0, 

which implies that r — tk+i — tk is the time between two consecutive captures. 

This paper follows a different approach. Indeed, we introduce a new sequence 
{tk}k of harvest instants, where the time for the [k + l)-th capture is determined 
by the amount of the biomass harvested at fc-th time (which will be denoted by 
x{tk)). This leads to a length of closed season determined by 

(1.2) tk+l-tk=T{x{tk)), 

where t(-) is some function depending on the amount of the biomass captured at 
t — tk. In consequence, the biomass harvested at time t = tk will determine the 
instant of the next capture. The idea is to define a length closure such that a bigger 
capture leads to a longest closure and by assuming broad and realistic properties 
on per capita growth rate and product function, we will find sufficient conditions 
ensuring a sustainable production, i.e., the existence of a globally stable periodic 
trajectory. 

This approach has been introduced in , where it is pointed out that the result- 
ing model is a new type of IDE equations, namely IDE with impulses dependent of 
time (IDE-IDT) and an introductory theory is presented. 

1.3. Outline. In section 2, we construct a model of fishery with closed seasons 
and pulse captures, which is described by an IDE-IDT. In section 3, we study 
some basic properties of the resulting IDE-IDT. The main results concerning the 
sustainability of the resource are presented in section 4. A numerical example is 
presented in section 5. 

2. The model 

A classical mathematical model of a fishery with closed seasons has to describe 
the following bioeconomic processes: i) The growth of the resource {e.g., a single 
marine species), ii) The production function, iii) The length of the closures. 
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2.1. Natural rate of growth. The growth of the biomass in the close season will 
be described by the ordinary differential equation: 

(2.1) x'{t) = x{t) r{x{t)), for any t € {tk, tk+i), 

where {tk} denotes an increasing sequence of harvest instants. 
Growth hypotheses (G) 
(Gl) Density— dependence. The per capita rate of growth r: [0,K] [0, +oo) 
is a derivable and strictly decreasing function of the biomass. In addition, 
we assume r(0) = tq > and r{K) = 0. 
(G2) Bounded variation. The rate r: [0,K] [0,+oo) has lowerly bounded 
derivative, i.e., there exists p > 0, such that ~p < r'(x) < 0. 

Remark 1. The property of derivability stated in (Gl) is a technical assumption. 
Nevertheless, in [25, some statistical results support the negative correlation be- 
tween per capita rate of growth and biomass. Notice that, and K are the unique 
equilibria of (|2.1|) and [0, K] is a positively invariant set. 

Remark 2. (G2) implies that the function x i— > xr{x) satisfies the local Lipschitz 
condition when x G [0, K]. Hence, the solutions of (|2.ip are unique to the right and 
depend continuously on initial condition to the right. 

An example of growth rates satisfying (G) are given by the family: 

(2.2) r(x)=ro(l- '^Y ^ ^^^^ ^^d l3 > 1, 

which generalizes the logistic Verhulst-Pearl per capita rate (see e.g., |29)V 
Another example is given by the function |26) : 

K ~ X 

(2.3) r{x)^ro— , with c> 0, 

K + cx 

which was used, for example, to describe the growth of Daphnia Magna. 

2.2. Fishing mortality. The fishing mortality (see [21] p. 102]) is the fraction 
F G [0, 1] of average population taken by fishing. In order to estimate it, we 
emphasize that there are two possible scenarios: an open season one, where the 
fishing is allowed all the time, and a restricted process (closed season), where the 
fishing is forbidden. 

The global fishing process is studied with two time scales: the first one describes 
the close season and only considers the growth of biomass summarized by (|2.ip with 
assumptions (G). The second scale concentrates the fishing mortality by considering 
the captures as pulses defined in a sequence of harvest instants {tk}, obtaining: 



(2.4) F 



x{tk)-xit+) _ H{x{tk)) 



x(tk) x{tk) 

where x{t'^) is the after fc-th capture biomass and F[{-) describes the capture func- 
tion in the biomass level at tk. From (|2.4|) . we can see that fishing mortality can 
be seen as a measure of catch per unit of biomass (CPUB). 

For convenience, let us define the impulse action /: [0,i4r] i— t- [0,/^], as follows: 

(2.5) I{x) = x~ H{x) ^ {1- F)x. 

In order to relate fishing mortality with the input factors (capital and labor) de- 
ployed along the fishing process, the continuous modeling literature has introduced 
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the concept of fishing effort, which is a rate describing the number of boats, traps, 
hooks, technicians, fishermens, etc., per time (see e.g., [1], [6], \21\). 

In an impulsive modehng framework, if the punctual fishing effort is denoted by 
i5 > 0, a bounded scalar measure, it is expected to describe (|2.4I) by a functional 
relation F — (f){E, x{tk)), which, combined with (|2.1I) and (|2.5I) . allows to introduce 
the complementary evolution law: 

(2.6) x{t+) = I{x{t)) = {l- (l){E,x{t)))x{t), with t = tk. 

We point out that the practical estimation of E and (j>{-, ■) are complicated 
matters in bioeconomic theory and we refer the reader to [T], [B], [30) for details. 
Harvest hypotheses (H): 

(HI) Impulsive action. /(•) is a derivable and increasing function. 
(H2) Elasticity. If Ax 0, then: 

, , AH/H H'(x)x 

for any x G (0, K]. This is, a percentage change in the resource biomass de- 
termines a percentage change bigger than or equal in the captured amount. 
Notice that, if the yield elasticity respect to the biomass is bigger or equal 
than one, then /(•) is inelastic or unitary. 

Remark 3. The property (HI) combined with (|2.6p says that: 
x^{E,x) <l-(j){E,x) for any x€[0,K]. 

Remark 4. The property (H2) says that a fixed punctual fishing effort is more 
productive at higher resource availability. In addition, by using (j2.5p , we can prove 
that (H2) is equivalent to 

I' ( x^ X 3d) 
\ \ < 1, and ■^(E,x)>0, for any xe(0,K]. 
I[x) ox 

Remark 5. Notice that, the called Cobb-Douglas production function can be in- 
terpreted by a fishing mor tality (j){E,x) = qE^'x^^''^ with g > 0, a > and ;9 > 0. 
It is easy to see that hypotheses (H) are verified when /3 > 1. An important 
case is the Schaefer assumption [25], corresponding to P — \ and a — 1, i.e., the 
parametrization is linear with respect to the fishing effort and biomass. 

2.3. Length of the closures. There exist a third evolution law governing the 
dynamics, which determine the sequence {tk\ of harvest instants. It is the first 
order recurrence that follows: 

(2.8) Atk=tk+i-tk=T{I{x{tk))), 

where r : [0, K*] [0, +oo), with K* ^ I{K). 

As it can be observed in (|2.8p . the length of the next closure, namely Atk, is a 
function r(-) of the stock after the fc-th harvest and allows to establish an automatic 
regulation of the dynamics by closed seasons. Here, we only introduce a theoretical 
proposal and we do not deal with the problems of implementation, for instance, 
those relating to the estimation of data requirred to define the length of the closed 
seasons. 
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2.4. The equation model. The dynamics determined by the combination of equa- 
tions (|2.1|) . (|2.4|) and ()2.8|) is formahzed by the impulsive differential equation: 

r x'{t) = x{t)r{x{t)), t^tk, 

(2.9) <^ xit+) = I{xit)), t^tk, 
[ Atk = r(/(x(tfc))), fc>0, 

where {t,x) G [0,+oo) x [0, +cx2). 

This type of impulsive differential equation is denoted as Impulsive Differential 
Equations with Impusive Dependent Times (IDE-IDT), which have been intro- 
duced by Cordova-Lepe in [9] and its novelty with respect to classic impulsive 
differential equations is that the sequence of impulse instants is determined by the 
process dynamics: the harvested stock I{x{tk)) will determine the next harvest 
time tk+i- Indeed, the sequence of harvest times is described by: 

(2.10) =<fc+T(/(x(tfc))), 

where the biomass x{t) is abruptly reduced to x{t^) = x{t) — I{x{t)) at t — tk- 

There exists several models of pulse harvesting of a renewable resource (not 
uniquely restricted to fisheries) described by impulsive equations, e.g.,: |^, g], [5], 
[S] and [32] consider a resource with logistic growth rate, [3S] considers a generalized 
logistic growth. In addition, Gompertz models (which, not satisfy (G2)) have been 
studied in [5], [M]. Nevertheless, these works consider a fixed time between two 
harvest processes, which is equivalent to consider r(-) as a constant function. 

We point out that impulsive models having sequences similar to ()2.10p have also 
been introduced in |20) by Karafyllis in an hybrid control theory setting and are 
named hybrid systems with sampling partition generated by the system. 

3. The impulsive system (|2.9p 

Given a first harvest time to G M and a biomass level xq e [0, iiT], then the 
existence, uniqueness and continuability of the solution of (|2.9|) . with initial condi- 
tion {tQ,xo), can be deduced from [9]. Indeed, we know any solution is a piecewise 
continuous function having first kind discontinuities at the harvest instants t = tk 
[k = 0, 1,2,---). In addition, we point out that different initial conditions will 
determine different sequences of harvest instants. 

3.1. Basic properties. In the study of (|2.9p . it will be necessary to consider the 
initial value ordinary associated problem: 

r z'it) = z{t)r{z{t)), 
^ > \z(<To)=w, with (cro,w) e K X [CiC]. 

Definition 1. The unique solution of p.l|) will be denoted by t t-^ Lp{t\aQ,v), for 
any t > ao, and ip{ao; ao, v) = v. 

Observe that given v G [Ojif], the function Lp{-; aQ,v) : [ao,+c>o) [0,K] satis- 
fies: 

(3.2) ip{t;ao,0) = and ^(t; (Jo, > 0. 

Let x{-) be the solution of (|2.9p with initial condition {to,xo), which determines 
the sequence {{tk,x{tk))}k- Since (|2.9I) is an ODE on {tk,tk+i], we can deduce that 
x(-) coincides with the unique solution (f{-,tk,I{x{tk))) of p.ip on (^^,^^+1]. 
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By using Definition [l] it follows that 

(p{a + tk;tk,I{x{tk))) = I{x{tk))expy J r[(p{s;tk, I{x{tk)))] ds 

= I{x{tk))exp J r[f{s + tk;tk,I{x{tk)))]dsy 

Finally, uniqueness of solutions implies (p{s + tk] tk, I{x{tk))) = ip{s; 0, 1{x{tk))), 
which leads to: 

(3.3) f{a + tk; tk,I{x{tk))) - I{x{tk)) exp ( / r[(p{s; 0, 1{x{tk)))] ds) , 



for any a G [0, r(/(a;(tfe))]. 

By using (|2.10p . it follows that at ct = T(I{x{tk))) {i.e., at t = tk+i), the solutions 
of (|2.9p satisfy the one dimensional map: 

(3.4) x{tk+i) = fix{tk)), 

where the function /: [0,K] — > [0,A'] is described as follows: 

(3.5) f{x)=F{I{x)), with i^(y) = y exp(^ *''V[(^(s;0,2/)]ds). 

Notice that Eq. (|2.6l) implies /(O) — 0. This fact will have important conse- 
quences when studying the asymptotic behavior of (|2.9p . 



3.2. Special solutions of ()2.9p and bio-economic interpretation. Let us in- 
troduce the straightforward result: 

Lemma 1. Any positive fixed point u* G (0,A'] of the map ()3.4p defines a piecewise- 
continuous T (I (u*)) -periodic solutiont u*{t) of (j2.9p . namely, the u* -associated 
solution. In addition, the fixed point u* — of the map (|3.4p defines a constant 
null solution of (|2.9p . 



The existence of a r(/(7i*))-periodic solution can be interpreted as a fishery 
strategy with harvest instants uniformly distributed in time. There exists different 
stability definitions for these solutions (see e.g., [9] and [20]). In this context, we 
will follow the ideas stated in [10| : 

Definition 2. The u* -associated solution of p.4p is locally asymptotically stable if 
there exists 6 > such that for any solution t ^ x{t) of (j2.9p with initial condition 
x(0) satisfying \x{Q) -- u* \ < 5, it follows that: 

lim \x{tk)-u*\^Q, 

k—y-\-oo 

where {tk} is the corresponding sequence of harvest instants associated to x{0). 



Definition 3. The u* -associated solution of (|3.4p is globally asymptotically stable 
if for any solution t i— >■ x{t) of ()2.9p with initial condition x{0) > 0, it follows that: 

lim \x{tk) -u*\ = 0. 

Observe that the asymptotic stability of a T(/(u*))-periodic solution implies the 
ecological sustainability of the fishery. On the other hand, the asymptotic stability 
of the null solution implies the future resource extinction. 
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4. SUSTAINABILITY CONDITIONS 

4.1. General result. 

Theorem 1. Let us assume that (G),(H) and the closed season hypotheses: 

(CI) Growth type. The function r: [0, i^*] — ^ [0, +oo) is derivable and de- 
creasing, such that t{K*) = 0. 
(C2) Initial condition. The initial value tq = t(0) satisfies 

(C3) Main condition. The following inequality: 

(4.2) |t(z2) - t(zi)| < 1 jln (^|^ - m{z, - z,) 

is verified for < zi < Z2 < I{K) = K* and m — p[e"'^° — 1]/q^; with 
a = pK + ro . 
are satisfied. 
Then: 

i) If €'^"'^"{1 — (/)(-£, 0)) < 1, then the resource-free solution of (12.91) is globally 
asymptotically stable [extinction case). 

ii) If e'^°'^"{\ — (f)(E,0)) > 1, then there exists a unique initial condition x* = 
f{x*) G {0,K) -with /(•) given by (j3.4p - defining a t (I (x*)) -periodic glob- 
ally asymptotically stable trajectory (sustainable case). 

Remark 6. Notice that (C3) gives us a range of graphic possibihties for the choice 
of function t(-). Moreover, (j4.ip impHes that the right side of (|4.2p is non negative 
for any z £ [0, K*]. 

Remark 7. If (C2) is verified, we can see that e'^"^°{l — (j){E,0)) > 1 imphes 
1| 1„(1 - ,(E. 0))| <r,<^^ 1„ (l + TT^(i + 5=)) ■ 

In consequence, the left inequahty says that there exists a trade-off between the 
fishing effort E and the maximal lenght of a closure Tq ensuring the fishery sus- 
taintability. In addition, observe that the inequality stated above has sense only 
when (l){E,0) G [0,4'*) C [0,1], which add a complementary constraint for the 
fishing effort. 



Proof. The asymptotic behavior of (|2.9|) is determined by the map p.4|) . 
Now, we will verify that / : [0, iiT] ^> M satisfies the following properties. 

a) The map / is derivable and /(O) = 0, 

b) The map / is increasing and f{K) < K , 

c) For any x S (0, K] it follows that: 

xf'ix) 

o< ^tV < 1- 

Indeed, a) is straightforward consequence from (|2.6p combined with f'{x) 
F'{I{x)) I'{x), where F'{y) is defined by: 

r / ^ I r / / \ ^ \i I / \ . I ^' r 



exp 



(/ r[^{s-0,y)]ds)[l + y(r[^{T{y)-Q,y)]T'{y)+ _ [^(s; 0, y)] dsj | 
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Let us verify that /'(O) = e^«''n/'(0) is consequence from (G1),(H1),(C2) and 
p.2p . Now, by using (HI), we observe that b) fohows if F'{y) > for any y e [0, K]. 
When dropping the exponential factor of F'{-), we only have to prove that 

(4.3) l + y(r[^(r(y);0,y)]T'(y)+y^ 0, y)] dsj > 
for any y G [0, K]. 

By hypotheses (G) and (CI), combined with f^(s; 0, y) > for any s € [0, t(?/)], 
we can observe that inequality (|4.3I) can be deduced from: 

(4.4) l>y(royiy)\+p 0, y) dsj . 

By integral representation of (p{s;0,y + h) and (p{s;0,y), with s e [0, t(2/)] and 
h > sufhcientl small, the Gronwall inequality implies that 

\^{s;0,y + h)-^{s;0,y)\ < \h\ e^^ ''+-<^> , 

for any s £ [0,T(y)]. Indeed, ||(s;0,y) < e(''^+'■o)^ with s e [0,T(y)]. Therefore, 
we can reduce our proof to demand the condition that follows: 

(4.5) y(£ e°'r{y)_^ -roT'(y))<l, where a^pK + ro- 

We point out that (C3) is equivalent to (|4.5I) . Indeed, when replacing zi and Z2 
by y and y + h {h > Q) respectively, (|4.5p is obtained by letting ft, -> 0. Inversely, 
since r(y) < tq = t(0), the inequality (|4.2p is obtained by integrating (|4.5p on 
[0,K% with X* = /(X). 

To prove that the right side of (|4.2p is greater than zero, we observe that 

Tn(ir*) - ln(z) 



inffi^^m^: zeiO,K*]} = ± 



and l/if* > m is equivalent to (C2). Finally, observe that (H) and (CI) imply 
f{K) = I{K) < K and property b) follows. 

The property c) is equivalent to x^{f{x)/x)' ~ F' {I{x))I' {x)x — F{I{x)) < for 
any x S (0, K]. This is verified if and only if for any x G (0, K], it follows that: 

(4.6) xl'{x) [1 + I{x)W'{I{x))] exp(T4^(/(a;))) < /(x) exp(W^(/(a;))) 

where W{I(x)) is defined by 



W{I{x))^ r[ipis;0,l{x))]ds, for xe[0,K]. 

Let us recall that 

rriy) g 

W'{y) = r[^{T{y);0,y)]T'{y)+ [^{s;0,y)] ds < 0, for ye{0,K*]. 

Jo 9y 

By using this inequality, combined with (j4.3p and Remark^ we can deduce that: 

(4.7) ^Trf <i< 1 ,LuT( ^i*^ ^e(0,i^]. 

l{x) 1 — I[x)\W'{l(x))\ 

By applying Lemma [5] (see Appendix), there are two posibilities for system p.4p 
according the sign of 

l-/'(0) l-e""'""/'(0). 
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By Eq. (|2.6p . we can express the threshold condition for case (a) /'(O) < 1 as 
e^or-o(i_0(£; 0)) < 1 and for case (b) /'(O) > 1 as e^" {I -(/){£, 0)) > 1. So that, 
the result follows. □ 

4.2. Application to logistic growth. Notice that, in some cases, the one-di- 
mensional map p.4p associated to the system ()2.9|) can be defined explicitly and 
the previous result improved. Indeed, let us consider a marine species with logistic 
growth, whose exploitation is described by: 

[ x'{t) = rx{t)[l-^), t^tk, 
(4-^) \ x{t+) - I{x{tj), t^tk, 

^ Atk = T{I{x{tk))), k>0. 

Corollary 1. Let us assume that the impulse action /(•) satisfies (H) and the 
close season satisfies (CI) and 

(C3') Closure condition. The following inequality: 

K..)-r(.,)l<^|ln(^)+ln(f^) 
ro I \zij \K - Z2/ 

is verified for < zi < Z2 < = K* . 



The 



//e''^"(l — 0)) < 1, then the resource-free solution of (|4.8p is globally 
asymptotically stable [extinction case). 

If eJ'^°{l — (f){E,Q)) > 1, then there exists a unique initial condition x* = 
f{x*) G (0,/^) defining a t [I [x*)) -periodic globally asymptotically stable 
trajectory {sustainable case). 



Proof. A simple computation shows that p.4[) is equivalent to the one-dimensional 
map: 

(«) .(.,.0 ^ /(..«,)) ^ ^ ,;,'!','gfo)')i,-,......,.», ■ 

We will verify that the map (|4.9p satisfies the assumptions of Lemma [2] (see 
Appendix). Firstly, notice that /(•) can be described as follows: 

f{x)^F{I{x)), with F{u)^——^\-—- 

u-\- [K — u\e '^^I'M 

and by using (H), it follows that /(•) is derivable and /(O) — 0. 

Secondly, observe that f'{x) — F'{I{x))I'{x). As in the proof of Theorem [1] 
it follows that F'{u) > (with < m < I{K) if*) if and only if (C3') is 
verified. By using this fact, combined with (HI), it follows that /(•) is increasing 
and f{K) = F{I{K)) < K. 

Finally, from (H2) combined with Remark SI it follows that 

g ^ xf^ ^ xI'{x)F'{I{x)) ^ xl'jx) F'{I{x))I{x) ^ F'{I{x))I{x) 



fix) F{Iix)) I{x) F{I{x)) - F{I{x)) ■ 

By using (CI), it is not difficult to show that uF'{u) < F{u) for any u G 
{0,I{K)]. This fact, combined with the last above inequality implies that < 
xf'{x)/f{x) < 1 for any x G (0,^]. Now, as /'(O) = e'^^''(l - (^(S, 0)), the result is 
a direct consequence from Lemma [5] □ 
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Remark 8. 

i) The assumption (C3') is equivalent to the differential inequality K > 
—ru{K — u)t'{u) for any < u < I{K), which furnishes a way to design 
admissible functions t(-) describing the lengh of open seasons. 

ii) In addition, it is easy to verify that in this case there are no explicit re- 
strictions for tq as stated in (C2) (see also Remark IH]). 

5. Example 

Let us consider a fishery with: biomass growth described by the logistic equation, 
fishing mortality satisfying Schaefer assumption, i.e., (j)(E,x) — qE and closures 
having lengths determined by the linear function r: [0, K*] — J> [0, +oo): 

(5.1) t{z) = a(K* - z), with a > 0. 

By using remarks [T] and [5l it follows that hypotheses (G) and (H) are satisfied. 
In addition, observe that assumption (CI) is satisfied since t(-) is strictly decreasing 
and t{K*) = 0. Moreover, notice that: 

—ru{K — u)t'{u) = aru{K ~ u) < arK^ /A, 

and by using statement i) from Remark HI it follows that (C3') is satisfied if 
a < A/{rK). 

On the other hand, as K* = (1 — qE)K and tq — r(0) = a(l — qE)K, we obtain 
the threshold: 

(5.2) £{qE) = e''^°(l - (j){E, 0)) = e'^''^^^^'?^) (1 - qE), 

and by Corollary [U it follows that the resource-free solution of (|4.8p is globally 
asymptotically stable \i £{qE) < 1. Similarly, there exists a globally asymptotically 
stable periodic solution if £{qE) > 1. 

In order to illustrate some properties of the set of punctual fishing efforts {0,E*) 
{£{qE*) = 1) ensuring sustainability, let us represent the slope of (|5.ip as follows: 

4 

a = rj, with < ry < 1. 

rK 

Notice that, to find E* is equivalent to find the unique fixed point w* — 1 — qE* 
of the map w t-^ e^^^'™. In this case, E* is dependent of the parameter ij G (0, 1), 
which is positively related to tq. 

It is straightforward to verify that the function rj i— E* is increasing and concave. 
This implies that lower values of maximal closure length tq leads to narrow ranges 
of sustainable punctual effort (0,i?*). 

We illustrate this previous results by using the numerical methods developed by 
Del-Valle [12]. The following parameters are employed: 

(5.3) ro = 0.05 [time^^], K = 1000 [tons], q = l and 77 0.25. 

These values determine w* 0.567143 with a respective E* - 0.432857. The 
figure shows the biomass curves associated to the following punctual fishing efforts: 



(5.4) 



Effort 
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Case 2 


Case 3 


Case 4 
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0.4329 


0.5 
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Biomass vs Time 
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Figure 3. Numerical solution of (|4.8p with Schaefer assumption, 
closure length defined by (|5.f |) and parameters (|5.3p . The punctual 
fishing efforts ensuring sustainability satisfy E E (0,0.432857). 

The Figure 3 shows the evolution of the biomass by considering an initial con- 
dition a;(0) = 300 tons and four different punctual fishing efforts. As stated before, 
we can see an extinction scenario for any punctual effort bigger than E* (this is 
the case for E = 0.5). Finally, it can be observed that the lenght of the closure is 
an increasing function of the punctual effort E. 



The mains results (Theorem [T] and Corollary [T]) propose sufficient conditions 
ensuring the ecological sustainability of a simple fishery model with impulsive cap- 
tures, which can be seen as the trade-off between the fishing effort E and the 
maximal closure length tq. The novelty (compared with harvest instants uniformly 
distributed in time) is to allow a variable length of closures, which has social and 
economic consequences in a short term. 

The proof of Theorem [T] is carried out by constructing a one-dimensional map 
(|3.4p . whose asymptotic behavior inherits crucial features of the IDE-IDT equa- 
tion. In this context, this approach could be extended in several ways. First, notice 
that assumptions (C2) and (C3) imply that the map (|3.4p is strictly increasing. 
Whe think that it is possible to consider more general maps and obtain less restric- 
tive conditions ensuring convergence towards a fixed point. This remains a future 
problem and its main difficulty is the hybrid nature of IDE-IDT equations. 

Provided that the fishery is sustainable, an important open problem is to find 
first order conditions on the punctual fishing effort E that maximizes: 



the sustainable production per unit time. In other words, the task is the Maximun 
Sustainable Yield (MSY) problem. 

Another extension of this work would be to consider the logistic case (|4.8p . 
replacing r > by a continuous function r : K — s> M+ . This problem is interesting 



6. Discussion 



(P{E,x*{E))x*{E) 



T{[l-q^{E,x*{E))]x*{E)) 
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by a bioeconomic point of view since periodic and Bohr almost periodic functions 
provide a good tool in order to modeling birth rates with seasonal behavior. 

On the other hand, to consider the logistic model (|4.8p with almost periodic 
perturbations has mathematical interest in itself. Indeed, if r(-) is a positive Bohr 
almost periodic function and r(-) is a constant function, it can be proved that 
there exist a unique almost periodic solution (in the sense of Samolienko and Per- 
estyuk ^24^), which is globally asymptotically stable. Nevertheless, there are neither 
equivalent results nor a study of asymptotic properties when r(-) is not constant. A 
careful development of the qualitative theory for IDE-IDT seems essential to cope 
with this kind of problems. 



Appendix 

The following result plays a key role in the proof of Theorem [1] 
Lemma 2. Let us consider the one dimensional map: 

(6.1) Xn+l = fiXn), Xo G [0, A'], 

where the function f: [0, AT] H> [0,/(Ar)] C [0, AT) satisfies the following properties: 

a) / is derivable and /(O) = 0, 

b) / is increasing and f{K) < K , 

c) For any x G (0, AT] it follows that: 

r. xf'(x) ^ 

o< < 1. 

fix) 

i) ///'(O) < 1 then it follows that lim Xn = 0. 

ii) // /'(O) > 1 then there exists a unique fixed point x* G (0, K) a it follows 
that lim x„ = x* for any Xq G (0, A']. 

Proof. Case /'(O) < 1: we can verify that there exists 5 > such that f{x) < x for 
any x G (0, 6). Let us denote by A, the set of positive fixed points of /. HA — 0, 
by using continuity of /(•) it follows that f{x) < x for any x G (0, AT]. 
If A 7^ 0, let us define ^ — mi{A}. It is straightforward to verify that: 

(6.2) ^>0 and f{x)<x for any x e {0,^). 

Now, by using c) we can deduce that < /'(^) < 1, which implies the existence 
of (5 > such that f{x) > x for any x G — 6,(_), obtaining a contradiction with 

In consequence is the unique fixed point in [0, K] and f{x) < x implies that 
{xn} is strictly descreasing and lowerly bounded. Finally, the convergence towards 
a; = follows from uniqueness of the fixed point. 

Case /'(O) > 1: We can verify that there exists S > such that f{x) > x for any 
X G (0, S). On the other hand, f{K) < K combined with continuity of / imply the 
existence of a fixed point x* G (0, AT) minimal with this property. Hence, it follows 
that f{x) > X for any x G (0,a;*) and (c) implies that < f'(x*) < 1. 

In order to prove the uniqueness of x* , let us define g{x) — f{x + x*) — x* and 
observe that g: [0, AT - x*] ^ [0,f{K) - x*] C [0, AT - x*) and .g'(0) < 1. Hence, 
the non-existence of fixed points on (x*,K) and the property f{x) < x for any 
X G {x* ,K] follows as in the previous case. 
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If xq G {x* , K] it follows that {xn} is strictly descreasing and lowerly bounded by 
X*. On the other hand, if xo G (0, x*) it follows that {x„} is strictly increasing and 
upperly bounded by x*. The convergence towards x* follows from the uniqueness 
of the fixed point in {0,K). □ 

Example: The function /: [0,K] — > R defined by f{x) — Xx/{a + x) with a > 1 
and X < 2a satisfies straightforwardly the properties a) and b) by choosing K > a. 
Finally, observe that: 

xf'(x) 1 

< i\' = < 1, for any x > 0, 

j(x) a + x 

and c) is verified. 

Hence, if A < a {i.e. /'(O) < 1), it follows that the sequence {xn}n defined 
recursively by (I6.ip with xq £ (0,iir] is convergent to 0. Otherwise, if A G (a, 2a) 
{i.e, /'(O) > 1), then the sequence is convergent to x* = a — A. 
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